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Abstract 

Theoretical investigations of the hyperfine structure of the hydrogen molecu- 
lar ion (one electron and two protons) are discussed. The nuclear spin-rotation 
interaction has been found to be of the same sign as in the hydrogen molecule 
and the hyperfine transition frequencies can be accurately predicted. With 
measurements of the hyperfine structure of the deuterium molecular ion (or 
of HD + ) it should be possible to obtain a value of the deuteron quadrupole 
moment that could be compared with the values obtained from the deuterium 
molecule and from nuclear theory. 
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INTRODUCTION 



The calculation of hyperfine structure (Ms) frequencies presents a substantial challenge 
to theory, in part because the experimental spectroscopic data available are so precise. For 
the hydrogen atom in the Is state, with electron spin S = | and nuclear spin / = |, the 
transition frequency between the F = and F = 1 states, where F = I + S, is known 
to be 1420.405 751 766 7(9) MHz Q, in comparison to the most complete theoretical result 
of 1 420.405 l(8)MHz 0, which includes reduced mass, QED, and hadronic (mainly proton 
recoil and size) effects. Indeed, a simple "Fermi" theory of the hfs (see, for example, |§, 
p. 74) with a fixed proton and the electron and proton ^-factors of, respectively, 2.002 319 
and 5.585 694, gives about 1422 MHz and with the correction due to the reduced mass 
of the proton-electron system it yields 1420.5 MHz. Further inclusion of the known QED 
corrections yields a value of about 1420.452 MHz Q and an additional ~ —45 kHz and 
the theoretical error comes from hadronic effects ||. Needless to say, the difficulty of the 
calculations increases substantially at each stage ||. 

The "hydrogen atom" of molecular physics is the hydrogen molecular ion H2 + consisting 
of two protons bound by a single electron. From a theoretical point of view this molecule 
offers a good starting point for investigations of molecular hfs since the electronic wave func- 
tion can be calculated to high accuracy. Surprisingly, there have been only a few experiments 
on the hfs of H2 + 0-|§, one experiment on HD + H^JT^] , and no experiments on D2 + . 

With the deuterated ions it should be possible to determine the electric quadrupole mo- 
ment of the deuteron, once the hfs transition energies are measured, and a further discussion 
of the prospects of using D 2 + to do this will be given below. 



THEORY OF H 2 + HFS 

As with the hydrogen atom, the hfs in H 2 + arises primarily from interactions between 
the electron spin and nuclear spin, but there are also smaller magnetic effects arising from 
the electron and nuclear spins interacting with the magnetic field of the rotating nucleus. 

The hfs energy can be described by the effective spin Hamiltonian 

H Ms = bl • S + d z S z + dS • N + fl ■ N, (1) 

where I, S, and N are, respectively, the total nuclear spin, electron spin, and rotational 
angular momenta, and b, c, d, and / are the coupling constants. Transition frequencies 
are obtained by diagonalizing the matrix elements of (jl]) with constants calculated from 
theoretical expressions or if given measured frequencies the constants can be determined 
by fitting using matrix elements of ([!]); usually the bps [0] angular momentum coupling 



scheme is used in which the intermediate vector F 2 = I + S is formed, and then the total 
F = F 2 + N. 

Early theoretical studies of H 2 + hfs were motivated by a need for predictions of transi- 
tion frequencies of the N — 1 rotational level of the v — vibrational state for astrophysical 
searches (it is now believed that there is little prospect for detection, cf. Hl2]), because no 
experimental data were available. The researchers calculated coupling constants and tran- 
sition frequencies using ([!]) with / = 0. Transition frequencies were estimated by Burke [13 
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using calculations of b, c, and d, and by Mizushima |14| who used values of b and c calculated 
by Stephens and Auffray [p~5H . A significant advance in the theory was the derivation by 
Dalgarno, Patterson, and Somerville |16|] of the b, c, and d terms of (|1|) from the Dirac eq 
for H 2 + by a nonrelativistic reduction within the Breit-Pauli approximation. They demon- 
strated that d could be written d = d\ + d 2 , where d\ is the electron spin-rotational magnetic 
interaction and d 2 is a second-order term in the electronic wave function. They calculated 
values for b, c, and d±, estimated d 2 , and calculated transition frequencies and the constants 
b and c were used by Richardson, Jefferts, and Dehmelt to make the first experimental esti- 
mate of d |J. Further improvements in calculations of b, c, and d\, and improved estimates 
of d 2 were presented by Somerville ]I7,I8| and by Luke ||19|| . The second-order constant 
d2 was calculated accurately using a variation-perturbation method by Kalaghan and Dal- 

who also calculated b, c, and d\ 



garno and by McEachran, Veenstra, and Cohen 
and accurate transition frequencies. 

The calculated constants discussed so far were all determined in the Born-Oppenheimer 
approximation in which the electron is assumed to move in the field generated by averag- 
ing over the electronic wave function. Nonadiabatic effects on the ground electronic state 
through the excitation of other electronic states due to the coupling of electronic and nuclear 
motion contribute significantly to the constant b and the effects of the reduced mass of the 
electron on the constants have been treated [^2,23]. 

The nuclear spin-rotation term / was introduced into the effective spin Hamiltonian ([!]) 
on phenomenological grounds by Jefferts [[/]], who measured the hfs transition frequencies 
for the v = 4-8 vibrational states. For the N = 1 rotational level for each v, he obtained 
values of the constants by fitting to ([I]). Varshalovich and Sannikov obtained different 
fits to Jefferts' data, notably, they found the value of / to be negative. The first (and to 
date the only) measurements of the hfs of the N = 1 level of v = state were carried out 
by Fu, Hessels, and Lundeen [[J through spectroscopy of high angular momentum Rydberg 
states of the hydrogen molecule H2 and by fitting the measured Rydberg transition levels to 
an effective spin Hamiltonian including ([!]) they determined values for b, c, d, and /. 

A recent theoretical study of the nuclear spin-rotation interaction in H 2 + found the sign 
of / to be negative |[24|| . The nuclear spin-rotation interaction can be written 



f(R) = f 1 (R) + f 2 (R), 



(2) 



and the major contributions [25|] to the energy hf are from the interaction of each nuclear 
magnetic moment with the magnetic field generated by the other rotating nucleus, 



hfx(R) 



R* ■ 



and with the magnetic field generated by the orbiting electron |26 



hf 2 (R) 



a- 



'R 2 



(3) 



(4) 



where in atomic units R is the internuclear distance, /i^ is the nuclear magneton, and the 
dimensionless quantities a, g p , and <7m are, respectively, the fs constant, the proton ^-factor 
defined previously, and the high-frequency component of the magnetic shielding constant as 
defined in Eq. (3) of |4|. 
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In Table [I] values of the hfs constants calculated as in p3| and the nuclear spin-rotation 
constants calculated as in [24|, averaged over the vibrational wave functions, are compared 

for the N = 1 level of the v = state and with constants 
on the N = 1 level of the vibrational state v = 4. 

-the theoretical constants in 



with constants measured 
determined [0 from the measurements 

The agreement between theory and experiment is impressive- 
Table | reproduce all measured transition frequencies to within 175 kHz. The formalism 
of |23[ accounts for reduced mass effects on the hfs constants but not fully for radiative and 



relativistic effects. 



DEUTERON MOLECULAR ION 

The electric quadrupole moment Q of the deuteron can serve as a sensitive test of models 
of the neutron-proton nuclear force f27|. The most accurate determination Q = 0.2860 ± 



0.0015 x 10~ 26 cm 2 was obtained semiempirically from measurements J^J] of the magnetic hfs 
of the deuterium molecule D2 through a theoretical value for the electric field gradient at the 
nucleus due to the other constituent charges of the molecule [ P5| , |50]] . It should be possible 
to measure Q semiempirically using D 2 + , which has the theoretical advantage of a simpler 
electronic structure than D 2 , thereby facilitating accurate computation, but to date there 
have been no measurements of the hfs of D 2 + . (An experiment that measured hfs levels of 
HD + was reported, but the quadrupole effects were not resolved |§.) 

For D2 + , the effective spin Hamiltonian is (|I|) plus an additional nuclear quadrupole 
interaction ]3TJ] proportional to eqQ, where e is the proton charge and q is the electric field 



gradient. Predictions of the hfs constants and electric quadrupole coupling constant for the 
N = 1 level of the v = state of D 2 + are given in Table ||. The constants b, c, and d 
have been calculated as in |23|] with the inclusion of nonadiabatic effects in b and c, and / 



has been calculated using the theory of p4| . The term eqQ was estimated by calculating 



q in the Born-Oppenheimer approximation and using the semi-empirical value of Q from 
D 2 . More detailed calculations of the hfs constants, quadrupole coupling constants, and 
transition frequencies will be reported Q . 

The quadrupole coupling constant is about 57 kHz, see Table [TI|. The experimental 
precision achieved in the ion trap experiments on H 2 + was better than 3 kHz (7|, in the 
Rydberg H 2 experiment it was around 20 kHz ||, and in a laser- rf double-resonance ion 
beam study of N 2 + it was better than 10 kHz It would appear that the prospects are 
good for measuring the deuteron quadrupole moment using D 2 + . 
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TABLES 



TABLE I. Comparison of EL" 1 " hfs constants from theory with those obtained from experimental 
measurements of hfs levels, for the N = 1 rotational levels of the v = and v = 4 vibrational 
states. All values are given in MHz and numbers in parenthesis represent errors quoted on the last 
digit. 

b c d f Ref. 



Theory, v = 880.163 128.482 42.421 -0.042 ]2|j24 

Experiment, v = 880.187(22) 128.259(26) 42.348(29) -0.003(15) g 



Theory, v = 4 804.104 98.008 32.658 -0.036 H|j24] 



Experiment, v = 4 804.087(2) 97.930(2) 32.649(2) -0.034(2) 111 



TABLE II. Predictions of the hfs constants and the electric quadrupole coupling constant for 
the molecular ion D2 + in the N = 1 rotational state of the v = vibrational level. All values are 
given in MHz. 



b 


c 


d 


/ 


eqQ 


135.739 


19.944 


21.461 


-0.002 


0.057 
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